Damping the zero-point energy of a harmonic oscillator 
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The thermal and zero-point energy of a damped harmonic oscillator are calculated for a range of damping 
values from zero to infinity. Both the thermal and zero-point energies decrease with damping. The energy 
stored in the quantum damped oscillator at fixed temperature increases with damping, an effect that may be 
experimentally observable. As the results follow from canonical quantization, the uncertainty principle is valid 
for all damping levels. 
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The one-dimensional damped harmonic oscillator has a 
central place in physics due to its simplicity and wide range 
of application. Although the classical treatment of this dy- 
namical system is elementary and its application to classical 
behaviour is well understood, the opposite is true in quantum 
physics. The dissipation of energy leads to difficulties in ap- 
plying the standard quantization rules to the damped oscilla- 
tor |T-5ri, and experimental damped quantum oscillators are 
probing the degree to which quantum mechanics describes the 
macroscopic world Il22lj26l . 

The technical difficulties in the canonical quantization of 
the damped harmonic oscillator can be overcome by inclu- 
sion of reservoir degrees of freedom that take up the dissi- 
pated energy. If the reservoir has a finite, or countably in- 
finite, number of degrees of freedom then a delicate limit- 
ing procedure must be employed to capture a wide range of 
damping behaviour |3, 10-20 1 (this limit is usually dealt with 
briskly, but the subtleties involved are lucidly demonstrated by 
Tatarskii 1 13 1). The limiting procedure amounts to a transition 
from a countably infinite to an uncountably infinite number 
of dynamical degrees of freedom in the reservoir In |21 J a 
reservoir with an uncountably infinite number of degrees of 
freedom 1271 was used from the outset, thus removing any 
need for a limiting operation. The resulting theory has many 
similarities with macroscopic electromagnetism, which can be 
viewed as a field theory of damped oscillators. In particu- 
lar, the dynamics of a general damped harmonic oscillator is 
governed by an arbitrary susceptibility that obeys Kramers- 
Kronig relations |21|,just as for the dynamics of light in a ma- 
terial medium. Canonical quantization of macroscopic elec- 
tromagnetism for arbitrary dielectrics |28 29 1 and canonical 
quantization of a general damped harmonic oscillator [21 J can 
both be directly achieved using the powerful reservoir formal- 
ism introduced in ll27l . 

In this paper we apply the general results for a quantum 
damped harmonic oscillator in thermal equilibrium |21] to a 
simple example and show that the zero-point energy of the 
oscillator is less than hu>Q/2, where ujq is the free oscillation 
frequency in the absence of damping. The energy removed in 
cooling the oscillator from temperature T > to its quantum 
ground state increases with damping, which offers one possi- 
bility of experimentally demonstrating the effect derived here. 

Following Ii2rj we consider a damped harmonic oscillator 



with unit mass, free oscillation frequency ojq and position op- 
erator g(t). The canonical momentum Ilq (t) is equal to the ve- 
locity and the canonical quantization relation [q{t),]Iq{t)] = 
ih is satisfied ||2TI . In thermal equilibrium the expectation 
values of the squares of the position and momentum operators 
are [21 1 
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where G{w) is a Green function 
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determined by an arbitrary susceptibility xi'^) that obeys 
Kramer-Kronig relations. When the damping is zero the sus- 
ceptibility x(w) vanishes and the Green function ^ is that 
of a free oscillator of frequency wq. The energy of the oscil- 
lator in thermal equilibrium must be obtained from the total 
thermal energy of the coupled oscillator/reservoir system. A 
simple prescription for performing this calculation was used 
in [21 1, and amounts to subtracting from the total thermal en- 
ergy the thermal energy of the reservoir in the absence of any 
coupling to the oscillator This gives the result 1,21 J 
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A similar procedure in the electromagnetic case gives the 
Casimir (zero-point plus thermal) stress-energy of the elec- 
tromagnetic field in a material 1281 . In the Appendix we show 
that the prescription giving (J4]| is equivalent to the thermal av- 
erage of the Hamiltonian of mean force [30|, consistent with 
earlier work on the thermodynamics of strongly coupled sys- 
tems fyi\. In light of this, the results in [28 1 show that the 
Casimir energy density is the thermal average of the Hamilto- 
nian of mean force for the electromagnetic field in a macro- 
scopic medium. 



The susceptibility x(w) of an experimental damped oscil- 
lator must be measured rather than postulated, just as the 
electromagnetic susceptibility of an individual material sam- 
ple must be measured (and will vary even for samples of the 
same material). The position and momentum correlation func- 
tions |21 1 of the damped oscillator in thermal equilibrium pro- 
vide one method of experimentally extracting x(w) and ojq 
that appear in Q. It is nevertheless useful to consider sim- 
ple formulae for the susceptibility of a damped oscillator, as 
is widely done in the electromagnetic case. In 121] the exam- 
ple of damping proportional to velocity was treated in detail, 
but this gives some problems with the zero-damping limit at 
T > if the corresponding susceptibility x{^) is taken to 
hold strictly at all frequencies up to infinity. The susceptibil- 
ity ceases to be physically meaningful at very high frequencies 
and therefore the high-frequency part of x(cj) should have no 
physical consequences. In addition, the case of damping pro- 
portional to velocity would not be expected to be experimen- 
tally relevant |21 ). We therefore here consider another exam- 
ple, chosen to be simple enough to allow exact analytical so- 
lution while being well-behaved in the limit of zero damping 
of the oscillator. 

We consider the susceptibility 
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where 71 and 72 are positive real constants. Being analytic on 
the upper-half complex-frequency plane, the real and imagi- 
nary parts of Q exhibit Kramer-Kronig relations and we can 

write 
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where a{Ld) is the function that couples the displacement q 
of the oscillator to the reservoir ETl . For (J5]l the coupling 
function is 
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The thermal results ([TJ, (|2]) and (|4]) were derived in 121] for 
cases where the total Hamiltonian is diagonalizable, a suffi- 



cient condition for which is the inequality 
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For the coupling function dTJi the condition dSll yields 

^0 > 27172, 
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so we impose this in what follows. When 72 — > 0, 71 re- 
maining fixed, the susceptibility (J5]) vanishes and the case of 
an undamped oscillator is recovered. 

The Green function pi for the susceptibility dSll has poles 



■'^-ni, W = -i72±wl, 
cji = Juj^ - 72(271 + 72). 
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The constants 71 and 72 thus serve as damping constants of 
the oscillator, while uji is a modified oscillation frequency 
when it is real. Note that (J9]) implies that the poles ( lOi are 
all in the lower-half complex-frequency plane so the Green 
function has retarded boundary conditions. The over-damped 
case with imaginary oJi can also occur while satisfying (|9]). 

The thermal expectation values (fTh, Q and (4|i can all be 
evaluated analytically for the susceptibility dSi. The inte- 
grands in each case are even functions of w for T > and 
so can be rewritten with lower integration limit of —00; the 
integrals are then evaluated by closing the integration contour 
in the upper (or lower) half-plane. The infinite sum over the 
residues of the poles of the hyperbolic cotangent function can 
be expressed in terms of harmonic numbers, but the resulting 
expressions are rather lengthy. In the limit 72 — ?> 0, the ex- 
pectation values Q, Q and Q reduce to the free-oscillator 
values {h/2ujQ)coth{hujo/2kBT), {hujQ/2)coth{hujo/2kBT) 
and {fkLjf)/2)coth{fkU()/2kBT), respectively (the momentum- 
squared expectation value is equal to the thermal energy for a 
free oscillator with unit mass, but for non-zero 72 these two 
quantities are different). The thermal energy as a function of 
the damping 72 is plotted for temperature T — fiu>o/kB in 
Fig.|2] 

The zero-point (T = 0) values of (fTl), (J2]) and (Bl) are most 
easily calculated directly rather than as the T ^- limits of 
the thermal results. The zero-point results are 
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We recall that these expressions presuppose the inequality (|9]l and note that they are real in the over-damped case where oJi 




FIG. 1 : Plots of the zero-point position uncertainty Aq (square root 
of l |12[ () and momentum uncertainty Ap (square root of (fl3l) versus 
damping 72 witii u)o = 10^" s^^ and 71 = ajo/4. The squares of 
the uncertainties are scaled with an appropriate power of Wq to have 
the same units as the product AqAp. The uncertainty relation is 
satisfied for all parameters obeying ^. 
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FIG. 2: Plots of the energy of the harmonic oscillator versus damping 
72, with LJo = 10"' s~^ and 71 = a;o/4, for T = and T = 
hujo/kB- The zero-point energy (T = 0) is damped below the free- 
oscillator value fujJo/2. The energy for T > is also damped below 
the free-oscillator value, though this damping is not very apparent 
except for very low T. The energy that can be extracted from the 
oscillator at T > (i.e. the T > energy minus the zero-point 
energy) increases with damping. 



is imaginary. 

Since all results here are derived from canonical quantiza- 
tion, the position/momentum uncertainty relation holds for the 
quantum damped oscillator. The squares Aq^ and Ap^ of the 
position and momentum uncertainties are given by ([12]) and 
( 13 I, respectively. Figure [T] shows the uncertainties as a func- 



tion of 72 for Wo = 10 s~ and 71 = ujq/A. As damping (72) 
increases the product Aq Ap increases from the minimum al- 
lowed value h/2. The position uncertainty Aq increases with 
damping while Ap decreases. 



the energy at temperature T — hujo/kB is also plotted. As 
72 increases the zero-point energy is damped below the free- 
oscillator value fiwo/2. In our example the oscillation fre- 
quency for 72 > is Wi, given by (111, provided cji is 
real. But the zero-point energy of the damped oscillator is 



not huii/2, as is clear from ( 14 1. In fact the range of 72 in 
Fig. l2]passes through wi = and into the over-damped case 
where ui is imaginary. The energy at T > is also damped 
below the free-oscillator value but the difference between the 
T > energy and the zero-point energy increases with damp- 
ing. This shows that the energy stored in the oscillator at fixed 
T increases with damping, and this energy must be removed 
to bring the oscillator to its ground state. 

We can also consider the limit of infinite damping, which 
occurs when 72 ^- 00. The condition (J9|l then requires 
71 -> 0, which we can satisfy by setting 71 = ajo/(472) for 
example. With this value for 71, the zero-point energy ([14]) 
goes to zero for infinite damping 72 — > 00, with a leading 
term of 
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The zero-point energy (14 1 of the oscillator is plotted in 
Fig.Sas a function of 72 for wq — IQ^^ s^^ and 71 = wo/4; 



The asymptotic approach of the zero-point energy to zero as 
72 —>^ 00 is thus very slow. We must note however that the 
position uncertainty diverges as Aq ^ 2\JhrfipK Jujq in this 
infinite-damping limit (with 71 = a;g/(472)). A large dis- 
placement of the oscillator can be expected to lead to nonlin- 
ear behaviour, so our assumption of a linear oscillator is not 
realistic for extremely large damping with susceptibility (]5]). 

As noted at the outset, the susceptibility is a quantity that 
must be measured, and in addition the "free-oscillation" fre- 
quency Wo is a parameter that must also be fitted to exper- 
imental data |!2T|. The theoretical ideal is a oscillator whose 
damping can be tuned from zero to a desired level, but this is a 
heavy demand in practice. A more realistic scenario is a set of 
macroscopic oscillators prepared with slightly different mate- 
rial geometries so that the damping varies slowly across the 
set. In the absence of data for the susceptibilities of such a set 
of oscillators, the results for the simple susceptibility (]5]l give 
some qualitative indications. The results in Fig. ]2] suggest that 
if all oscillators in the set are brought to a fixed temperature, 
then the energy removed from the oscillators in reaching their 
ground states will increase with damping, where the damping 
level is determined from the measured susceptibilities. 

We are indebted to J. Anders for informing us of the Hamil- 
tonian of mean force. We also thank N. Kiesel and M. As- 
pelmeyer for interesting discussions. 



Appendix 

Consider a system composed of two interacting parts; the 
system of interest (5), and a reservoir (i?). The total Hamilto- 
nian of this system is of the form, iJ = Hs + iii + Hr, where 
Hi characterises the coupling (of arbitrary strength) between 



S and the reservoir. We ask the question, what is the energy 
of S in thermal equilibriuml 

A choice of Hamiltonian, H* that gives the correct equihb- 
rium properties for S without reference to R is the Hamiho- 
nian of mean force 13^ 
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where ^ = l/kT, and Zr = Tri^[exp (-/3iJfl)]. The parti- 
tion function Z* associated with H* is then 
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where Z — Tr[exp (— /3iJ)] is the total partition function. It is 
evident that equilibrium averages of quantities pertaining to S 
alone, computed using H* will be identical to those calculated 
from the full Hamiltonian H. The factor of Z~^ within the 
logarithm plays no role in such a calculation, but is determined 
by the requirements that (a) when iJ/ -^ then H* — > Hs; 
and (b) the free energy F* associated with S is llSTl 



^* = -r^ iog(z*) = f-fr 
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which is the amount of energy available to do work in a re- 
versible, isothermal change of state of S, including that ob- 
tained through decoupling S and i? 131]. In the case consid- 
ered in the main text, where the g-oscillator plays the role of 



S, (A.18 I will give the correct generalized force (and there- 
fore work done during any isothermal change of state) when 
F* is differentiated with respect to the "free-oscillation" fre- 
quency wq^ or the quantities 71 2 within the coupling of the 
oscillator to the reservoir. Furthermore, when 71 2 — > then 
F* ^ Fs. 



Given the above properties, (A. 16 1 is interpreted as the ef- 
fective Hamiltonian of S in thermal equilibrium. In answer to 
our initial question, the equilibrium average of the associated 
energy is, using (A.17[), 
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where (FIr) = —d \og{ZR)/dl3 is the equilibrium average of 
the energy of R in the absence of any coupling to S. This is the 
prescription that was previously used to calculate the Casimir 
energy density [28 1 and the thermal energy of a damped har- 
monic oscillator II21I . the latter of which is given by (Hh. 
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